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A variety of methods are developed for characterising the free energy landscapes of continuum, Landau-type 
free energy models. Using morphologies of lipid vesicles and a multistable liquid crystal device as examples, 
I show that the methods allow systematic study of not only the most relevant minimum energy configurations, 
but also the transition pathways between any two minima, as well as their corresponding energy barriers and 
transition state configurations. A global view of the free energy landscapes can then be visualized using either 
a disconnectivity graph or a network representation. Different forms of free energy functionals and boundary 
conditions can be readily implemented, thus allowing these tools to be utilised for a broad range of problems. 
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I. INTRODUCTION 

In many soft matter systems the stability of minimum free 
energy configurations and the transition pathways from one 
state to another play a central role. Liquid crystal devices of¬ 
ten exploit the existence of two or more (meta) stable states 
given the device geometry and surface anchoring |[i|-[l], and 
are of great technological interests because power is only 
needed when the devices switch states. A major obstacle in 
many applications of superhydrophobic surfaces is the lack 
of stability of the suspended states nil, where gas cush¬ 
ions are trapped in between the surface corrugations. Shape 
transformations of lipid bilayers hold the key to understanding 
important issues in cellular transport (e.g. budding phenom¬ 
ena, encapsulation of nanoparticles) as well as how cells react 
to external stimuli ffl. All these examples require knowl¬ 
edge on the free energy landscapes, yet such calculations are 
usually limited to a certain facet of the landscapes, idealised 
scenarios, or both. Furthermore, while established methods 
exist to compute free energy barriers and transition states for 
atomistic simulations iflOl - U^ . equivalent methods at the con¬ 
tinuum level HEl are still lacking. 

It is therefore my aim here to describe a suite of methods 
capable of fully characterising the free energy landscapes of 
continuum models, including basin-hopping global optimiza¬ 
tion ilillii, doubly-nudged elastic band lO, and a hybrid 
eigenvector-following 111^ techniques. These methods will 
allow us to calculate not only the minimum energy states, but 
also the transition state configurations, the height of the en¬ 
ergy barriers, and the transition pathways between local min¬ 
ima. In some cases, the analysis will in fact show that there 
could be several competing pathways. It is also worth not¬ 
ing that, while these methods are well-known for atomistic 
systems, this paper takes the important step of reformulating 
them in a suitable form for continuum systems. 

I will focus on cases where the physical systems of interest 
can be described using Landau theory. The free energy of the 
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system is thus written in terms of expansion in scalar, vector 
and/or tensor fields. To highlight the strength and versatil¬ 
ity of the methods, I will discuss examples from the field of 
soft matter: (i) shape transformations of lipid vesicles and (ii) 
multi-stable nematic liquid crystal devices. It is straightfor¬ 
ward to implement other Landau free energy models which 
are routinely used in other areas of physics, chemistry and 
materials science. Complex boundary conditions (including 
Neumann, Dirichlet, periodic) can also be accounted while 
keeping the core of the simulation methodologies unchanged. 

This paper is organized as follows. The details of the com¬ 
putational methodologies and implementation are described 
in sections II and III. The methods are then used to study two 
different problems in soft matter and the results are presented 
in section IV. I summarise the most salient points of this paper 
in section V. 


II. THE FREE ENERGY FUNCTIONALS 

Conceptually the most important step in the implemen¬ 
tation of energy landscape methods is to define the energy 
(scalar) function to be optimized in terms of the degrees 
of freedom. In particle-based simulations, the degrees of 
freedom are usually the positions and orientations of the 
molecules. In field theories, they are the values of the fields at 
given coordinates in space. To illustrate how the free energy 
functionals can be defined in a suitable form for the methods 
described in the section III, I will now consider two specific 
examples. Here I shah assume it is appropriate to discretize 
space in a uniform grid, and for simplicity, I will choose a 
square grid in two dimensions and a cubic grid in three di¬ 
mensions. 


A. Phase Field Model for Lipid Vesicles 

To model the shape of lipid vesicles, I employ the phase 
field formulation proposed by Du and co-workers In 
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this model, the total free energy is given by 

^ = ^b + ^kv{V-Vof + ^kA{A-Aof. (1) 

The first term is the bending energy, 

2(02 _ 1)0^ dv, (2) 

where a scalar order parameter 0(x) is used to distinguish the 
inside (0=1) and outside (0 = —1) of the vesicle, e is the 
width of the interface, and k is the bending rigidity. The sec¬ 
ond and third terms in Eq. o are penalty terms to conserve 
the vesicle volume and area. Vq and Aq are the target volume 
and area, and ky and kA are constants. I use the following 
definitions to compute the volume and area of the vesicle. 


The above free energy functional must be discretized and 
written in terms of the relevant degrees of freedom. The goal 
is to describe the free energy and its derivatives as functions 
of {0ij/c}, where the subscript i{= 1... Nx), j{= 1... Ny), 
and k{= 1... Nz) label the lattice points in three dimensions. 
The bending energy, vesicle volume and area can be expressed 
as 
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and similarly for the derivatives in the y and z directions. 
Another possibility is to approximate dx(l)ijk = (0(i+i)j/c — 
0(2 -i)j7c)/ 2. However, this stencil is prone to a numeri¬ 
cal checkerboard instability if the bending energy term is re¬ 
moved. All the the stencils used here are second order accu¬ 
rate. 

Landau-de Gennes Free Energy - The liquid crystal appli¬ 
cation to be considered in this paper is a device first proposed 
by Tsakonas et al. in. The device consists of an array of 
square wells filled with nematic liquid crystals, and the well 
surfaces are treated such that the liquid crystal molecules lie 
parallel to the plane of the surfaces. This system is effec¬ 
tively two dimensional. Using the Landau-de Gennes frame¬ 
work iH, the liquid crystal configuration can thus be mod¬ 
eled by a symmetric traceless tensor. 


Q = 


Qll Qi2 
Qi 2 —Qll 


s{n (g) n — I), 


( 8 ) 


where n is the director and s is a scalar order parameter that 
measures the degree of orientational ordering. 

In the absence of external fields and surface effects, the 
Landau-de Gennes free energy under the one-constant approx¬ 
imation is given by [HI] 


Tr = 




-aTrQ2 + :^(TrQ2)2 


dA. 


(9) 

The first integral corresponds to the elastic energy of the liq¬ 
uid crystals, and tZei > 0 is an elastic constant. The second 
integral is the bulk free energy. The parameter c is material- 
dependent constant. Given T* is the nematic-isotropic transi¬ 
tion temperature, I define the coefficient a = 7 (T* — T), with 
7 > 0 and T is the temperature of the system. Here I shall 
assume T < T* such that a > 0. 

(5) To reduce the number of input parameters in the model, I 
will now rewrite the free energy functional in its dimension¬ 
less form. Substituting Eq. ©to (0 and defining x = x/L, 
^ kei = KeijalA, and ja, I ob¬ 


tain 


+ AxAyAz. 

ijk 

Specific stencils are needed to approximate the derivatives. 

Lor the Laplacian, I use 

(A0)yfe = {dl(l))ijk + {dl4>)i3k + {dl4>)i3k, (6) 


with 


{dx^fyijk — 
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In this form, there is only one input parameter in the model, 
the dimensionless elastic constant hei • Note that L is the size 
of the system, and from now on, I will drop the tilde in the 
equations. 

Using a square lattice, the discretized total free energy may 
be written as 

4 / = [\VQu\l + \VQ32\l]AxAy 


and similarly for the y and 2 ; directions. Lor the square- 
gradient term, one possibility is given by 


+ ^(Q?i+Q?2-l7Aa:A2/. (11) 


{dx4’fi3k 


[(0(»+l)jfc ~ 4>ijk)^ + (0(i-i)jfc — (t>ijk)‘^] 

2(Aa;)2 


( 7 ) 


7' 

Since the form of the gradient terms in Eq. cni) is the same as 
in (|4]), the stencil described in Eq. (|7]) may be used. Thus, for 
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either Qii and Q 12 , 


B. Transition State Search 


IV7^|2 _ [(Q(i+l)i Qij) + {Q(i-l)j Qij) ] 

' “ 2(Ax)2 

- QijY + (Qi(i-i) - QijY] 

2 ( A ^)2 


Suitable boundary conditions are needed at the solid 
walls. Here I simply consider the strong anchor¬ 
ing limit, where the director n lies parallel to the 
walls. In particular, I will follow the formulation given 
in Luo et al. H where (Qii(a:,y),Qi 2 (x,y))wau = 
s(a;,y)(cos20(a;,y), sin26»(a;,y)), with 0(0, y) = 0(1, y) = 
7 r/ 2 , 9{x, 0) = 9{x, 1) = 0, and s{t, 0) = s{t, 1) = s(0, t) = 
s(l, t) = f {t). I define 


fit) = < 


tjd, 

1 , 

(1 -t)/d 


0<t<d, 

d<t<l-d, 

l-d<t<l, 


(13) 


with d chosen to be The strong anchoring boundary 

conditions are implemented as Dirichlet boundary conditions 
on the Q tensor, which can be applied via additional virtual 
lattice nodes along the walls. They influence the optimization 
of the free energy functional through the gradient terms. 


III. ENERGY LANDSCAPE METHODOLOGIES 
A. Einding Minimum Eree Energy Configurations 

One of the primary goals in energy landscape methods is 
to characterise potential/free energy landscapes in terms of 
the minima and transition states (or more generally, saddle 
points). 

To obtain relevant local minimum configurations (low-lying 
energy states), I apply the basin-hopping algorithm ifLA fl^. 
In basin-hopping, a step consists of a trial move followed by 
an energy minimization. The simplest trial move consists 
of random perturbations for the lattice field values, = 
4>ijk + and similarly for the Q-tensor. Here ^ is a random 
number between -1.0 and 1.0, and d is the amplitude of the 
perturbation, which I usually take to be d = 0.125.. .0.5. 
The energy minimization is carried out using the limited- 
memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) algo¬ 
rithm EH [mi. All that is required for the LBFGS algorithm 
is the gradient of the free energy with respect to the degrees 
of freedom, e.g. d'^/d(j)ijk for the phase field model of lipid 
vesicles. 

The move is then accepted or rejected based upon the 
change in the corresponding free energy for the local mini¬ 
mum, A simple approach, which has been widely used in 
the literature, is to use a Metropolis acceptance criterion. I fur¬ 
ther note that, for continuum systems, the number of minima 
is often very small (< 10). In this case, it could be preferable 
to accept all steps, and as such, search for the majority, if not 
all, of the possible minima. 


Methods for finding transition states (saddle points in the 
free energy landscapes) can generally be divided into single- 
and double-ended search methods. Here I first use a doubly- 
nudged elastic band (DNEB) method lIT^ to find candidate(s) 
for transition state(s) between every pair of minima. Note that 
any two pair of minima need not be connected by a single 
transition state. These candidates are then refined using a hy¬ 
brid eigenvector-following technique E3, which is a single- 
ended search technique. An alternative transition state finding 
method which has begun to find applications in continuum 
models is the string method 

Once the transition states are found, small displacements 
are applied in the two downhill directions. Energy minimiza¬ 
tions are carried out using the LBFGS algorithm MM. 
Such a run results in the minimum energy pathway for a 
minimum-transition state-minimum triplet. In some cases, 
new minimum energy states may be identified after this step. 
By iterating the process, I can systematically build a database 
of minima and transition states, and how they are intercon¬ 
nected. 

Doubly-Nudged Elastic Band - The first step in the DNEB 
approach is to make an initial guess for a set of images, 
between two endpoints (minima) and 
p(N+i) j symbol as a shorthand for all the de¬ 

grees of freedom in an image a. For the vesicle application, it 
corresponds to {(j)ijk}^, whereas for the liquid crsytal applica¬ 
tion, it refers to {{Qii,Qi 2 )ijk}^- I typically use 25 images 
between the two endpoints. 

The simplest way to obtain an initial guess is to use lin¬ 
ear interpolation of the lattice field distributions between the 
two minima (end points). While this is suitable for the liq¬ 
uid crystal application, my results indicate that such linear 
interpolation is inefficient for phase fields models. To illus¬ 
trate this point, consider a liquid droplet (or a vesicle). Here 
the scalar field (p tends to two distinct bulk values, +1 in¬ 
side the droplet and —1 outside, with a smooth transition at 
the interface over a finite width of several lattice spacings. If 
one linearly interpolates the phase field distributions between 
two shapes of droplets, there will be a large number of lat¬ 
tice points in the intermediate images, which are not at the in¬ 
terface, yet their phase field values deviate significantly from 
+1 and —1. These configurations are not physically relevant. 
Instead, one needs an algorithm that linearly interpolates the 
interfacial profile, and thus the droplet shape. 

One possible algorithm is as follows. Generally, any lat¬ 
tice point in the simulation box falls into one of three cat¬ 
egories. In the first category, the lattice point is inside the 
droplet/vesicle in both minima. I assign 4>ijk = 1 for all the 
images. In the second category, the lattice point is always out¬ 
side the droplet/vesicle and (j)ijk = — 1 can be assigned to 
all images. Further, the boundary profiles for these two cate¬ 
gories may be computed. They are represented by two series 
of points, {ri}in and {rijout- In the third category, the lattice 
point is inside the droplet/vesicle in one minima and outside in 
the other minima. To assign a suitable value of cjfijk in image 
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a, I compute the weighted distance squared 

^ijk) ^ijk^in’) 

<^?j/c,out = (1 + ^?jk)‘^^ijk,out ^ (14) 

where 0“.^^ = (1 - j^)4>ijk + ^HkM and 

dijk,out are respectively the minimum distance between the 
corresponding lattice point and the two profiles, {rijin and 

{ri}out. If dj^k,m < 4fe.out> f set = 1. Otherwise, 
(jyfjk = —1. At this point, every image is a binary map 
of 0 = zbl. One can further smooth the image by taking 
advantage the fact that the scalar field follows a hyperbolic 
tangent profile near the interface ia. 

The images are relaxed using energy gradients that include 
contributions from two components. The ''true'' gradient, 
depends on the derivatives of the the free energy functional 
with respect to the lattice degrees of freedom. The spring 
force component, keeps the images roughly equidistant 
and is derived from the spring potential 

= \ ((«“’■)' - («“’^)') • (15) 

and are respectively the distances to the a — 1 
(left) and a + 1 (right) images of image a, (s^’^)^ = 

I pail _ pa 1^ phe coefficient k may be taken as a constant. 
Alternatively, it can be dynamically adjusted based on the de¬ 
viation of the image spacing from uniformity. 

Using the true and spring gradients without modification 
leads to corner-cutting (images are pulled away from the min¬ 
imum energy path) and sliding-down problems (images slide 
down from barrier regions) [Til]. Certain projections for 
the gradients are needed based on the unit tangent along 
the path. I follow the formulation in lITlIl . where is de¬ 
fined based on the tangents to the left and to the right image, 
^a,± _ pa±i _ pa image a is at a maximum or a 

minimum. 


r = 


ra,+AUro 
^ ^ n 

U/ypro 


-A4'" 

^ rr 


where 


if ^a+i ^ 
if ^a+l ^ ’ 

(16) 


Removing the tangent part of the true gradient is numerically 
not sufficient to achieve equidistant images. An additional 
component from the spring force is required, and in DNEB, 
the following component of the spring constant gradient is re¬ 
tained lIT^ . 


Sdneb — ^ 

+ g" -(g" •gi)g", (19) 

where gJ — (g°^ • 

Eigenvector-Following Technique - To refine the transition 
states obtained in the doubly-nudged elastic band procedure, 
I use an eigenvector-following technique In 

this work, rather than computing the complete Hessian (sec¬ 
ond derivatives of the free energy functional with respect to 
the lattice degrees of freedom), I use a Rayleigh-Ritz ap¬ 
proach ll^ to estimate the smallest eigenvalue and the cor¬ 
responding eigenvector. Since it only requires the gradient 
information [[l9l|^, it is more efficient for large systems. 

It is also essential that the zero eigenmodes are removed. 
This can be achieved using projection operators |29|] or 
eigenvalue shifting d. In the models described here, the 
zero eigenmodes correspond to translations in lattice direc¬ 
tions where periodic boundary conditions are applied. For a 
scalar field, and assuming periodic boundary condi¬ 

tion in the a-direction, the zero Hessian eigenvector is given 
by 


J_ / 94)111 d4>ijk d4>N^NyN^ \ 
N \ da da da y 


( 20 ) 


where A" is a normalisation constant, such that |e| = 1. Pe¬ 
riodic boundary conditions are not used for the liquid crystal 
application considered here. In cases where periodic bound¬ 
ary conditions are needed, the zero Hessian eigenvector can 
be obtained as in Eq. (l20l) . with the scalar field 0 replaced by 
the Q-tensor. 


IV. RESULTS 
A. Vesicle Morphologies 


If the image a is neither at a maximum nor a minimum, then 
the tangent vector is defined as follows 

if ^a+l > ^a ^ ^a-1 

(17) 

if ^a+1 ^ ^ ^a-1 

The unit tangent vector is defined as 

The first projection in DNEB retains only the components 
of the true gradient that are perpendicular to the unit tangent 
vector , 

g" =g"-(g"-T")fA (18) 


In this sub-section, I will consider the morphological tran¬ 
sitions of single-component lipid vesicles. I shall assume that 
the spontaneous curvature of the membrane and the Gaussian 
bending modulus are zero. Thus, the shapes of the vesicle de¬ 
pend only on the reduced volume [H^], k' = ’ where 

Vo and Aq are the volume and area of the vesicles, u is tuned 
in the simulations by varying Aq, while keeping the volume 
Vo constant. Here I have used al40xl40x210 grid and 
the following simulation parameters: n = 1.0, e = 1.5Ax = 
0.03, kA = kv = 10^ and Vb = 1-375 x 10^Ax^. I will 
also restrict the analysis to axisymmetric vesicles, though the 
methods themselves are not limited by this symmetry. 

While the shapes of lipid vesicles have been analysed in 
much detail in the literature (see e.g. iS HEIM El), these 
studies are usually limited to minimum energy configurations. 
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FIG. 1: Shape transition from a prolate (b) to an oblate (d) for a 
vesicle with reduced volume v — 0.987. Panel (a) shows the energy 
profile for the minimum energy path. The pear-shaped configuration 
in panel (c) is the corresponding transition state. 


My aim here is to show that, using the methods developed 
here, it is possible to systematically study not only the min¬ 
ima, but also the transition state configurations and the mini¬ 
mum energy paths between various vesicle morphologies. To 
illustrate this, I first consider a vesicle with a reduced vol¬ 
ume of = 0.987. When v is close to 1, there are only two 
possible minimum energy shapes, a prolate (Fig. [Hb)) and 
an oblate (Fig. Wid)). Furthermore, I find that there is only 
one transition state, the morphology of which is given in Fig. 
[TJc). Thus, the free energy landscape is very simple. Fig. 
[TJa) shows the minimum energy path for the transition from 
a prolate to an oblate; it corresponds to the transition path¬ 
way with the smallest possible energy barrier. The x-axis is 
a measure of distance in the configurational space between 
morphology a along the path to a reference minimum con- 
figuration, (s")^ = \ A movie 

showing the full transition pathway is available in the support¬ 
ing information ll^ . 

When the reduced volume differs significantly from 1, the 
free energy landscape becomes more complex. For u = 
0.661, there are (at least) ten minimum energy shapes. These 
minima are obtained from a basin-hopping run of 500 steps. 
Solutions arising due to the periodicity of the grid are ignored. 
Doubly-nudged elastic band id and hybrid eigenvector¬ 
following methods |[T^ are then used to compute the transi¬ 
tion states and minimum energy paths for transitions between 
these various morphologies. 

With an increasing number of minima and transition states, 
one needs a better way to visualize the free energy land¬ 
scape. Plotting the energy profile as in Fig. [TJa) provides 
detailed information on the transition between any two min¬ 
ima. However, the ''global view'' of the landscape is not cap¬ 
tured. One possible representation is using the disconnectivity 
graph iHI], as shown in Fig. [2l Here the ^-axis corresponds to 
the energy scale. The end of every line corresponds to a min¬ 
imum. If one traces two minima, they meet at a vertex. This 
vertex corresponds to the highest energy along the minimum 
energy path for the transition between the two minimum con¬ 
figurations. Thus, a disconnectivity graph provides informa¬ 
tion on the competing minima and how easy for these minima 
to interconvert. 

I note that the minimum energy shapes shown in Fig. [2] 
have different genera {"number of holes"). Thus, in gen¬ 
eral, the Gaussian curvature term should be taken into account 



FIG. 2: A disconnectivity graph for a vesicle with reduced volume 
a = 0.661. 


when the transition pathways are computed. Extensions of 
the model described in section IIA to include Gaussian cur¬ 
vature, multicomponent membrane and adhesion to solid sub¬ 
strate are available in the literature (see e.g. imHHH), and 
can be readily implemented in the energy landscape methods 
described here. These are avenues for future research and are 
outside the scope of this report. 


B. Multistable Nematic Liquid Crystal Devices 

The second application to be considered is a multistable ne¬ 
matic liquid crystal device nil . The device consists of an ar¬ 
ray of square wells filled with nematic liquid crystals, and the 
well surfaces are treated such that the liquid crystal molecules 
lie parallel to the place of the surfaces. One such cell is mod¬ 
eled in Fig. fusing the framework of the Landau-de Gennes 
free energy, assuming strong anchoring on the surface. I use 
a 150 X 150 grid, and the dimensionless elastic coefficient is 
taken to be Kqi = d x 10“^. 

Within such a well, there are two classes of minimum con¬ 
figurations, namely diagonal and rotated states ||1,[1,[3^. Due 
to the symmetry of the system, there are two equivalent di¬ 
agonal states, see Fig. Oa-b), and four equivalent rotated 
states, see Fig. 13c -f). Since there are six minimum free en¬ 
ergy states, in principle there are 15 transition pathways to be 
considered. However, transitions between any rotated and di¬ 
agonal states are related by symmetry. Further, I find that the 
minimum energy path between any two rotated states or two 
diagonal states is composed by a sequence of rotated-diagonal 
transitions. A direct pathway between pairs of diagonal states 
or pairs of rotated states is never observed. Thus, the free en¬ 
ergy landscape of the system may be represented by a network 
shown in Fig. [S] Each line corresponds to a rotated state- 
transition state-diagonal state triplet. Such a network repre- 
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FIG. 3: Minimum energy configurations in a square liquid crystal 
device. There are two classes of configurations: (a-b) diagonal and 
(c-f) rotated states. Each double arrow corresponds to a minimum- 
transition state-minimum triplet, thus visualising the connectivity of 
the free energy landscape. 

sentation is an alternative to the disconnectivity graph used in 
the previous sub-section to visualise the landscape globally. 

For a given rotated-diagonal pair, there are in fact two com¬ 
peting transition mechanisms. As shown in Fig. lU^a), the en¬ 
ergy profiles for the two mechanisms are virtually the same. 
In the first mechanism, Fig. Ub-d), a bend defect transforms 
to a splay defect in the top left corner. This results in an excess 
— 1/2 defect which propagates to the bottom left corner, and 
subsequently accommodates a splay to bend defect transfor¬ 
mation there. In the second mechanism, Fig. HJe-g), the tran¬ 
sition is driven by a +1/2 defect, which propagates from the 
bottom left to the top left corner. The full pathways for these 
two mechanisms are available as movies in the supporting in¬ 
formation . It is worth mentioning that the transition states 
with defects shown here are reminiscent of recent experimen¬ 
tal observations in shallow nematic chambers S- Higher en¬ 
ergy pathways are also possible and are observed nimerically, 
corresponding to transition mechanisms mediated by two or 
more defects. This will be discussed in a separate publication. 

V. DISCUSSIONS 

Here I have described a set of methods suitable for map¬ 
ping the free energy landscapes of continuum, Landau-type 


free energy models in terms of their minima, transition states, 
and minimum energy paths. Compared to studies which focus 
only on the minimum energy configurations, this approach 
provides several significant advantages. It allows us to sys¬ 
tematically study the energy barriers and transition pathways 
as a function of the system parameters. On a fundamental 
level, it enriches any study which involves the construction 
of a phase or morphology diagram. As shown above, not 
only there are competing minima, but also competing path¬ 
ways may exist for the transition between two minimum con¬ 
figurations. It is also possible to analyze the (dis)appearance 
of minima, and correspondingly, how the connectivity of the 
landscapes is altered as a result of changing the system pa¬ 
rameters. On a more practical level, these energy landscapes 
techniques have great potential to be used as a design or op¬ 
timization tool. Taking the multistable liquid crystal devices 
as an example, the minimum energy path corresponds to the 
most efficient transition pathway, including the minimum en¬ 
ergy required to induce the transition between any two min¬ 
imum states. It can also provide insights into the effective 
form of external perturbations (e.g. strength and direction of 
the external electric field) required for the transition. Addi¬ 
tionally, since the configurations of the most relevant minima 
and transition states are known, these results may be exploited 
to design suitable surface treatments (anchoring strength, sur¬ 
face corrugations, etc), which increase/decrease the stability 
of the minima and alter the heights of the energy barrier for 
the transitions. 


It is worth noting that it is straightforward to implement 
other forms of free energy functionals. Given the widespread 
use of Landau free energies in physics, chemistry, and mate¬ 
rials science, these methods may therefore have wide-ranging 
applications. The descriptions and discretizations of the the 
free energy functionals are also compatible with other sim¬ 
ulation techniques, including the implementation of complex 
boundary conditions, such as surface corrugations, curved and 
non-symmetric boundaries. I envisage the use of these energy 
landscape techniques to be highly complementary to methods 
which capture the dynamics of the system (such as finite ele¬ 
ment, lattice Boltzmann method, etc). 
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